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1 Introduction 

The CDF and D0 experiments have successfully collected data since start of the Run 
II at the Tevatron Collider in 2001. The large B-meson production cross-section and 
the possibility to produce all kind of B hadron states, opened to the two collaborations 
the possibility to study with high precision the tiny effects of CP-violation in the 
Heavy B hadrons systém, and to search for new physics effects in rare decays, in a 
way unavailable to the previous generation experiments. A new and largely unknown 
sector of the Heavy Flavor physics, complementary to the one already tested with 
precision at the B-factories, as recently pointed out by I. Bigi p], has begun to be 
explored in search of possible signs of new physics. In this short note a selection of 
the most recent results on heavy B hadrons (mostly B s mesons) from the Fermilab 
Tevatron, and from the Belle experiment running at the Y(5>S) are reviewed. 

2 CP Violation in B s mesons at Tevatron 

In the neutral B s systém, the CP asymmetry in — > J/ip<ft decay play the analogous 
role of the B° — > J/ipK® for the B^ systém. Decays of the B s meson via b — > ččs 
transitions in fact can be used to probe, via interference effects in the mixing and 
decay amplitudes of the process, the f3 s angle of the squashed (bs) unitarity triangle, 
defined as (3 S = arg (— Vt s V t l/V cs V* b ). An important difference with respect to the 
Bd systém is that in the Standard Model (3 S is expected to be very small (~ 0.02), 
making the measurement of CP asymmetry in B® — > J/ifxfr decay very sensitive to 
new possible physics effects in the mixing phase of the B® — B s systém. 

CDF and D0 experiments have been able, for the first time, to perform a search 
for CP violation in the neutral B s meson systém, by measuring the time-dependent 
CP asymmetry in the B® — > J/ip4> decay mode. The vector-vector finál statě J/ip4> 
contains mixtures of polarization amplitudes: the CP-odd A±, and the CP-even A and 
A\\ amplitudes. These terms need to be disentangled, using the angular distributions, 
in order to extract /3 S , and their interference provides additional sensitivity [2]. 
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Figuře 1: Data samples for the CDF and D0 analyses of B® — > Jifxp decay. The CDF 
experiment reports 3166 ± 56 events, while the D0 experiment reports 1967 ± 65 
events. 



Very recently the CDF collaboration presented a new updated flavor-tagged, time- 
dependent, analysis, based on 2.8 fb -1 of integrated luminosity [3], that supersedes 
results from the previous 1.7 fb^ 1 untagged analysis [I], and 1.35 fb -1 flavor-tagged, 
time-dependent analysis [5]. The D0 collaboration result is also based on a 2.8 fb^ 1 
sample of flavor-tagged data, and is published in [6]. The B® — > J/ip4> signals from 
the two experiments are shown in Fig. [TJ 
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Figuře 2: Confidence regions in the space of parameters Ar s and (3 S for the CDF 
(left), and D0 (right) analyses. The green band corresponds to new physics models, 
as described in the text. 

In the CDF analysis two different fits have been performed, one without using 
flavor tagging and assuming the Standard Model, which fixes (3 S — 0, and a second 
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one with flavor-tagging and letting (3 S to be free. Assuming no CP violation (fl a = 0) 
the fit allows to simultaneously measure the decay width difference Ar s and the 
average lifetime of the B s meson. The results are reported in Table [U 



Table 1: Standard Model fits in the CDF tagged analysis. 



Parameter 


CDF measurement (untagged) 


CT S = 2c/(T H + T L ) 


459 ± I2(stat) ± 3(sysť) fjm 


AT S = T L - T H 


0.02 ± 0.05(síaí) ± 0M(syst) ps" 1 




A \ 2 


0.508 ±0.024 ±0.008 




^||(0)| 2 


0.241 ±0.019 ±0.007 



It is worth noticing here that the average B s lifetime measurement is consistent 
with the HQET expectation of equal lifetimes for the B s and Bd mesons. Table [1] 
reports also the measured polarization amplitudes for the B® — ► J/ip(p decay, that 
are consistent with those measured in the similar B^ decay B° — ► J/ipK*° [7|. For 
the CP fit, CDF does not report point estimates for any of the physics parameters, 
providing instead the confidence region in the {(3 Sl AT S ) plane shown in Fig. [2] (left), 
computed from Monte Carlo using Feldman-Cousins method for confidence intervals. 
Also shown in the same figuře, is the theoretical expectation from the Standard Model 
(black point), and in presence of new physics (green band). The Standard Model p- 
value calculated using the likelihood ratio is of 7%, corresponding to 1.8 Gaussian 
standard deviations. 

Treating AT S as a nuisance parameter, CDF reports confidence intervals for /3 S , 
and find that (3 S is within [0.28, 1.29] at 68% C.L., and within [-vr/2, -1.45] U [-1.01, 
-0.57] U [-0.13, vr/2] at 95% CL. 

To remove the two-fold ambiguity in the likelihood of the time-dependent, flavor 
tagged analysis, D0 constrained the strong phases of the helicity amplitudes in the 
B® — > J / ip(j) decay to the world average values for the B° — > J / ipK*° decay, measured 
at the B-factories [8]. Some justification of the constraint used has been recently 
showed in [3]. 

Results of the constrained fits are shown Table [2J In this case three types of fit 
have been performed: a Standard Model fit which fixes f3 s to its expected value, a CP 
fit with f3 s floating, and in addition, a CP fit with the further constraint that Ar s = 
2 1 Tf 2 1 cos (p s , where <p s = arg(—M[ 2 /ri 2 ), is the mixing phase of the B s syst and 
M( 2 is the off-diagonal element of the mass matrix governing the flavor oscillations 
in the B s systém (related to the mixing frequency by Am s = 2|Mf 2 |)- The average 

the discussion the approximation <f) s ~ — 2(3 S has been made. This is a reasonable approx- 
imation since, although the equality does not hold in the Standard Model, both are much smallcr 
than the current experimental resolution, whereas new physics contributions add a phase 4>NP to (f> s 
and subtract the same phase from 2(3 Sl so that the approximation remains valid. 
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Table 2: CDF results from tagged B° s —>■ J/ip(p analysis. 



Parameter 


CP Fit 

((p s floating) 


SM Fit 

(<l>s = 0) 


NP Fit 

(Ar s = 2 \M( 2 \ cos (p s constraint) 


n (ps) 


1.52 ±0.06 


1.53 ±0.06 


1.49 ±0.05 


Ar s (ps- 1 ) 


0.19 ±0.07 


0.14 ±0.07 


0.083 ±0.018 


A ± (0) 


0.41 ±0.04 


0.44 ±0.04 


0.45 ±0.03 


|4>(0)| 2 -|4l(0)| 2 


0.34 ±0.05 


0.35 ±0.04 


0.33 ±0.04 


(ps = -20. 


u - ' -0.30 


fixed (-0.04) 


-0.46 ±0.28 



lifetime r s , and the decay amplitudes are consistent with expectations and with the 
CDF measurements. Confidence regions in the ((p s (= — 2/3 s ), Ar s ) space are shown 
in Fig. [2] (right). The point estimate for the CP violation phase and Ar s , obtained 
by D0 are: (p s = -0.57íg;^(síaí)íg;^(sysí) and Ar s = 0.19 ± 0.07 (stat)^° f(sy st) 
ps -1 . The Standard Model p-value is of 6.6%, and the fTuctuation respect to the 
Standard Model goes in the same direction as CDF. 

Experimental sensitivity to new physics effects on the B s mixing phase can also 
be obtained from charge asymmetry measurements in semi-leptonic B s decays. The 
semi-leptonic asymmetry is in fact related to the mixing phase by the relation: 



AS N(B S -> f) - N(B S -> f) Ar 

A%r = — 7= : ; =r ~ tan<p s , 

SL N(B s ^f) + N(B s ^f) Am s 



where / corresponds to direct B s decays B s — > / (e.g. D~l + vi). 
The Standard Model prediction for the semileptonic asymmetry in B s decays is 
very small, at the level of few units in 10 -5 |10j . 

At the Tevatron the semileptonic asymmetry has been measured both in inclu- 

N + + —N 

sive di-muon samples, where A S SL ~ ^ M N >* M , or by using the sequential decays 

sample B° s — > ^ivD s - The first method has very high statistical accuracy, but requires 
knowledge of asymmetries of other contributing processes in addition to the detector 
charge asymmetries. The second method has less statistical power but ensures that 
the major contribution to the asymmetry comes from the B s decays. Combining the 
two measurements D0 obtains: A' SL = 0.0001 ± 0.0090 [lij, while the CDF result 
based on 1.6 fb -1 di-muon pairs, is: A S SL = 0.020 ± 0.028 [12J. At the current level of 
precision A S SL is not able to provide powerful constraints on new physics contributions 
on the mixing phase. 

Both the CDF and D0 analyses of the CP violation in B° s — > J/ip(p decay show a 
slight disagreement with the Standard Model prediction, and both results ffuctuate in 
the same direction. Recently the D0 collaboration has made public the results of the 
fit without the strong phase constraints, allowing the HFAG group to combine them 
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with the CDF results [13] . The combined countours are shown in Fig. [3j The p-value 
for the combined result is 3.1%, corresponding to 2.2 Gaussian standard deviations. 
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Figuře 3: HFAG Combinations of the CDF and D0 B® — > J/ip<p tagged analyses. 
The plot on the right use additional experimental input from the CDF and D0 
measurements of A% T . 



3 Rare B s decays at Tevatron and Belle 
3.1 Leptonic Two-body Decays at Tevatron 

The FCNC process — > fi + fi~ is predicted to have a branching ratio of B(B® — > 
= (3.42 ± 0.54) x 10" 9 in the Standard Model [T4], well below the current ex- 
perimental sensitivity of the Tevatron experiments. The B° — > /í + /í~ decay is further 
suppressed by |^d/V^ s | 2 , with a predicted branching ratio of (1.00 ± 0.14) x 10~ 10 . 
Significant enhancements are instead predicted by several new physics models. For 
example in the minimal super-symmetric standard model (MSSM) the B® branching 
ratio is proportional to tan 6 /3 where tan/3 is the ratio between the vacuum expec- 
tation values of the two neutral Higgs fields. In iž-parity violating super-symmetric 
(SUSY) models an enhancement is possible even at low values of tan/3. 

Both Tevatron experiments have dedicated triggers to collect B — > fi + fi~ events, 
and optimized selections for the B® — > fi + \i~ candidates based on sophisticated mul- 
tivariate analysis techniques. D0 combines the discriminant variables in a likeli- 
hood ratio, while CDF uses a neural network (NN) discriminant. Both experiments 
estimate the dominant combinatorial background by a fit to the mass sidebands, 
while contribution from decays of B mesons to two light hadrons, which could peak 
in the signál mass region, is estimated to be an order of magnitude lower than 
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the combinatorial background. Both experiments do not report any significant ex- 
cess of signál candidates over the expected background, and set 90% CL. limits 
calculated with a Bayesian method for a data sample of 2 fb _1 per experiment: 
B(B° S -> fi+fi-) < 7.5 x 10" 8 (D0 ) [15J and B(B° S -> n+fT) < 4.7 x 1(T 8 (CDF) [16]. 
Because of the superior mass resolution of the tracking systém CDF is able to sepa- 
rate B° s and B° mesons and to quote a 90% CL. limit separately for the -B decay of 
B(B° -> y+ir) < 1.5 x 10- 8 . 

Very recently the CDF experiment reported the results of a search for the B° sd ^ 
e + e~ channel and for the lepton-flavor violating mode B° sd — > e + fi~. In particular 
the latter decay modes are strongly suppressed within the Standard Model, in which 
leptons do not change flavor. These decays are allowed, however, in many models of 
new physics, such as Pati-Salam leptoquarks model, or in SUSY and Extra Dimension 
models, where the assumption of a local gauge symmetry between quarks and leptons 
at the lepton-flavor violation tree-level couplings leads to the prediction of a new force 
of Nature which mediates transitions between quarks and leptons [T7] . 

Using a 2 fb _1 data sample, CDF find one event in the search window for the B® — > 
e + fi~, with estimated 0.81 ± 0.63 background events, and one event for B° s — > e + e~ , 
with 2.8 ± 1.8 estimated background events. By using the -B — > K + n~ decay mode 
as a relative normalization, CDF derives the upper limits at 90% CL. on the decay 
branching ratios of B(B° S -> e+//~) < 2.0xl0~ 7 , and B(B° S -> e+e") < 2.8 xlO" 7 pj]. 
Finally from the decay branching ratio limits CDF calculate the corresponding lower 
bound on the Pati-Salam leptoquark mass: M L q(B° s ) > 47.7 TeV/c 2 at 90% CL. (see 
Fig.HD. 
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Figuře 4: Leptoquark mass limit corresponding to the 90% CL. on B(B® — > e + fi 



3.2 Charmless Two-body Decays at Tevatron 

Non-leptonic decays of b hadrons into pairs of charmless charged hadrons are effective 
probes of the CKM matrix, and sensitive to potential new physics efřects. The large 
production cross section of b hadrons at Tevatron, and the ability of CDF to trigger 
on fully hadronic decays, allows extending such measurements to and A° decays, 
complementing the B° meson case, extensively studied at the B-factories. 

CDF analyzed an integrated luminosity of ~1 fb^ 1 of pairs of oppositely-charged 
particles, selected by the displaced track trigger. A sample of 14500 H® — > h + h~ 
decay modes (where H® = B°,B® or A° and h = K,7i,p) was reconstructed after 
the off-line confirmation of trigger requirements. The invariant mass resolution and 
the particle identification separation power available in CDF, are not sufficient to 
disentangle the individual H® — > h + h~ decay modes on an event-by-event basis, 
therefore a Maximum Likelihood fit is performed to separate the different components. 
The fit combines kinematic and particle identification information, to statistically 
determine both the contribution of each mode, and the relative contributions to the 
CP asymmetries. 

Significant signals are seen for B° — > 7r + 7r~, B° — > K + n~ , and B® — * K + K~ , 
previously observed by CDF [19]. In addition to that, three new rare decay modes 
have been observed for the first time B° s — > K~n + , A° — ► pn~ and A° — > pK~, with 
a significance respectively of 8.2a, Q.Oa and 11.5c. No evidence was obtained for 
B° s -> 7T+7T- or B° -> K+K~ mode. 

Table 3: Branching fractions results. Absolute branching fractions are normalized 
to the the world-average values B(B° -> K + n-) = (19.4 ± 0.6) x 10" 6 and f s /f d = 
0.276 ± 0.034 and f A /f d = 0.230 ± 0.052 [18]. 



Mode 


BR(10- 6 ) 


B^ 


7T + 7T 


5.02 ± 0.33 ± 0.35 


B° s ^ 


K+K~ 


24.4 ± 1.4 ± 3.5 




K-7T + 


5.0 ± 0.7 ± 0.8 






5.6 ± 0.8 ± 1.5 






3.5 ± 0.6 ± 0.9 






0.49 ± 0.28 ± 0.36 


B ^ 


K+K- 


0.39 ± 0.16 ± 0.12 



The absolute branching fractions obtained by CDF normalizing the measurements 
to the the world average value B(B° — > K + tt~) = (19.4 ± 0.6) x 10~ 6 , are listed in 
Table [31 while the CP-related measurements are listed in Table IH where fd, f s and f\ 
indicate the production fractions respectively of B° B° s and A° from fragmentation 
of a b quark in pp collisions. 
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Table 4: CP-violation related results. 



Mode 


Measurement 


A C p(B° -> K+Ti-) 
A CP (B° S - K-n+) 

A CP {K ~+pK~) 
A CP {A° b ^pn-) 


-0.086 ± 0.023 ± 0.009 
0.39 ± 0.15 ± 0.08 
-0.37 ± 0.17 ± 0.03 
-0.03 ± 0.17 ± 0.05 



The decays A° — > pn~ and A° — > pK~ are allowed at tree level in the Standard 
Model, but are suppressed by the small value of the involved CKM matrix element 
V u b- Loop diagram processes can contribute at a magnitude that is comparable to the 
tree diagram process, leading to sizable direct CP violation. In the Standard Model 
a Aqp value of 0(10%) is predicted. 

The measurement of the direct CP violation asymmetries in the 6-baryon decays, 
presented by CDF, is the first such measurements in this sector. The statistical 
uncertainty still dominates the resolution and prevents a statement on the presence 
of asymmetry, whose measured value deviates from at 2. la level in the A° — > pK~ 
decay mode and is fully consistent with in the A° — > pir~ decay mode. In Fig. [5] the 
invariant mass spectrum in the tt + tt~ mass hypothesis, and the relative probability 
density function for the A° — > pK~ are shown, illustrating the good description of 
the data by the fit and the powerful A°/A° separation. 
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Figuře 5: Invariant mass spectrum for it + tt mass assignment (left) and relative 
probability density function (pdf) of A£ -> pK~: pdf(Ag)/[pdf(Ag) +pdf(Á~£)] (right). 
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3.3 Radiative B s penguins at Belle 

During the last several years the possibility of performing B° s meson studies at the 
e + e~ colliders running at the T(5S) resonance has been extensively explored. The 
first evidence for production at the T(5S) was found by the CLEO collaboration 
[26], [27] using a data sample of 0.42 fb -1 collected in 2003. This study indicated that 
practical B® measurements at the T(5S) are possible with at least 20 fb -1 , which can 
be easily collected at B-factories running with ~10 34 cm _2 sec _1 luminosity. To test 
the feasibility of a B® physics program the Belle collaboration collected at the T(5S) 
a sample of 1.86 fb -1 of data in 2005, and a one of 21.7fb _1 in 2006 [2Š1I29]. 

The collected samples have been used by Belle to perform searches for exclusive 
radiative decays of the B® mesons, that are of great interest, because sensitive to new 
physics effects and experimentally unaccessible to the Tevatron experiments ,(the 
presence of low energy photons in the finál statě makes these kind of decays too hard 
to be reconstructed in CDF and D0 ). In particular Belle searched for the decay 
modes B° s — > 07 and B® — > 77, using the full 23.6 fb -1 available data sample |30j . 

w+ 



Figuře 6: Diagrams describing the dominant SM processes for the B° s — > 07 (left) 
and Bg — > 77 (right) decays. 

Within the Standard Model (SM) the B® — > 07 decay is described by the radiative 
penguin diagram shown in Fig. [6] (left). The branching fraction is predicted to be 
~ 4 x 10 -5 |31J. The B® — > 77 decay proceed via the penguin annihilation diagram 
shown in Fig. [6] (right) and is expected to have a much smaller branching ratio, in 
the range (0.5 — 1.0) x 10~ 6 , that can be however enhanced by about an order of 
magnitude in various new physics models [321 (33j [34], reaching a level not far from 
the current sensitivity of the Belle experiment. 

To extract the signál yields a multi-dimensional un-binned extended maximum 
likelihood fit is performed to the M^ c and AE variables. Mb c and AE are respectively 
the beam-energy-constrained mass, and the energy diff erence observable, defined as: 
AE = E%f - E™ m and Mbc = v / (ECM m )2 _ ( p c^)2 j where E ™ and p CM are the 

energy and momentům of the B° s candidate in the e + e~ center-of-mass (CM) systém, 
and E^ m is the CM beam energy. 

Fig. [7J shows the M^ c and AE projections of the fit results of the data, together 
with the fitted functions. A clear signál of I8Í5 events is seen in the B® — > 07 mode, 
with a significance of 5.5 standard deviations, providing the first observation of a B s 
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Figuře 7: M^ c projection (left) and AE projection (right) for the B° s — > 07 (top) and 
— * 77 (bottom) modes. The thick solid curves are the fit functions (thin solid 
curves: signál functions, dashed curves: continuum contribution). 

penguin radiative decay. The branching fraction is measured to be B{B® — > 07) = 
(5.7tJ:f(stat.)í};?(syst.)) x 10~ 5 , in agreement with the SM predictions. 

No significant signál is observed instead for the B® — > 77 mode, and an upper 
limit at the 90% CL. of B(B° S -> <jry) < 8.7 x ÍO" 6 is set. 
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